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Abstract 

We present a parabolic approximation that incorporates reflection. With this approximation, there 
is no need to solve the parabolic equation for a coupled pair of solutions consisting of the incident and 
reflected waves. Rather, this approximation uses a synthetic wave whose spectral components manifest 
the incident and reflected waves. 
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The (Leontovich-Fock) parabolic approximation, which approximates the elliptic Helmholtz equation by 
a parabolic partial differential equation, was originally applied to electromagnetic propagation. EJ Tappert 
and Hardin introduced the parabolic approximation to acoustic propagation in the ocean to account for 
inseparable range effects in the sound speed profilc.0 In ocean acoustics, the parabolic equation is a useful 
computational tool for tackling inseparable indices of refraction for which the sound speed profile changes 
slowly with respect to range. 

One of the deficiencies of the standard parabolic approximation is that it neglects backscatter. Heretofore, 
to account for backscatter, one solved the parabolic equation for a coupled pair of solutions (the incident 
and reflected solutions). Attempts to account for backscatter include among others the works of Collins 
et al, which uses a two-way parabolic approximation.Ll Herein, we present a different approach to include 
backscatter. Based on rigorous ray tracing, we combine the incident and reflected waves into a modulated 
synthetic wave that progresses in the incident direction. 

Rigorous ray tracing has been developed in a generalized Hamilton- Jacobi representation that accounts 
for terms ignored by classical ray tracing and other asymptotic methods. It has provided insight into 
propagation phenomena. Rigorous ray tracing has shown that the existence of a sound-speed gradient 
is sufficient to induce linear (material) dispersion and angular (geometric) dispersion even for isotropic 
frequency-independent sound-speed profiles, that rays are not generally orthogonal to wave fronts, that 
classical ray tracing does not predict all caustics, and that rigorous ray tracing may be solved in closed form 
whenever the corresponding wave equation may be solved in closed form.Q Its quantum mechanical analogy, 
the trajectory representation, has shown how to derive the Helmholtz equation (the stationary Schrodinger 
equation) from the generalized Hamilton- Jacobi equation.l3 This allows us to construct the wave function or 
normal mode from Hamilton's characteristic function (a generator of the motion for the trajectory or ray 
path). These normal pcnodes can be synthetic normal modes that contain the incident and reflected waves as 
spectral components.El We shall use such a normal mode to develop the parabolic equation that accounts for 
reflection. 

Our objective in this letter is to present a parabolic equation that accounts for reflection. It is beyond 
the scope of this letter to solve the resultant parabolic equation. The acoustical community is free to solve 
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this equation by the methods of their choice. This work is presented in two dimensions, which is sufficient 
to illustrate how to incorporate reflection into the parabolic equation. 

We assume that the ocean to first order is a stratified medium whose index of refraction varies with depth 
due to temperature and pressure changes. The range dependence of the index of refraction is second order. 
This index of refraction is dependent upon two cartesian coordinates: (x, z) for range and depth respectively. 
The index of refraction varies much more rapidly in the z-direction then in the x-direction. We also assume 
that, for propagation of a wave train through the ocean medium, the reflected wave is much smaller than 
the incident wave consistent with the concept of backscatter. 

Recently, the trajectory representation of quantum mechanics (the quantum analogue to rigorous ray 
tracing) showed how the reflected and incident waves can be combined to synthesize a wave whose front 
monotonically advances in the direction of incidenceO. The synthetic wave is given by 
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wave front moves in +a;-direction 



where a is the amplitude of incident wave and /3 is the amplitude of reflected wave, where \P\ < \a\, and 
where k is the wavenumber and w is the angular frequency. This synthetic wave is a normal mode (follows 
from the superposition principle). 

The synthetic wave has spatially modulation in phase and amplitude as shown by right side of Eq. (|^) . 
For completeness, the right side of Eq. (^ was derived from the generator of the motion for the trajectory 
in Ref. and the left side was subsequently developed by analysis. While the right side of Eq. (|l|)|-.was 
first developed from Hamilton's characteristic function by the quantum analogy to rigorous ray tracing,Q we 
subsequently learned how to do it in a wave representation. This is the contribution of rigorous ray tracing 
that we use here. 

The wave equation in two dimensions is given by 

d^-^/dx^ + d'^-^/dz^ = C-^ix, z)dH/di^ 

The speed of sound, C, is isotropic and only spatially dependent. The wave equation is separable in time 
so that ^(x, z, t) = Tpix, z) exp(i[j;i) Hence, the wave equation is reduced to the two-dimensional Hclmholtz 
equation 



d^ij/dx^ + d^ij/dz^ + K^{x, z)ip = 

where k{x, z) — lo/C{x^ z). 

For reference, the standard parabolic approximation substitutes 
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'0(X, z) = 6'standard(a;)0standard(a;,2) = exp(ifca;)0standard (s^, 2) 

into Eq. (|^) to produce, after a standard simplification, the standard parabolic equation given byi 



dz"^ ' dx 
which does not incorporate reflection. 

Let us incorporate reflection by considering 



'standard , C'tpstandard , / 2 7 2\ ; n 
»2fc h (k - fc )0standard = 0, 
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tpix, z) = 9{x)(f>{x, z) 



where 
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0=[a^ +f3^ + 2aPcos{2kx)]^^^ exp 



,arctan(^tan(fc.x^ 

a + p 



There is flexibihty in choosing the form of 9. Different choices of 9 lead to different parabolic equations.i 
We have chosen a 9 that is the spatial component of the synthetic wave, Eq. (|l|). This 9 includes reflection 
while progressing in the incident direction. In the standard parabolic equation, the corresponding ^standard 
in Eq. ^ is given by ^standard = exp(jfca:;), which only includes the incident wave. Substituting — 9(1) into 
the Helmholtz equation leads to 

where d^9jdx^ = k^9 by the superposition principle or by direct substitution. 
We now examine {d9/dx)/9, which is given by 



89 /dx 
9 

For small reflections, /3 <C a, Eq 
d9/dx 



f3cos{2kx) 
13 cos(2fca;) 



1/2 



X exp i arctan ^^3^ tan(fca;)^ — i arctan ^^q^ tan(fca;) 
. may be simplified to 

= ik[l - {213 /a) cos{2kx)] exp[i(2/3/a) sin(2fcj;)] + 0[{[3laf] 
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Now the transformed Helmholtz equation for small reflection becomes 



ttI + ttI + «2fc[l - {2(3 /a) cos{2kx)] exp[i(2/3/a) sin(2fca;)] 
az ox 
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The critical assumption for the validity of the parabolic assumptions is that is well behaved (smooth) in 
range so that 



d^(j)/dx' 



2kd(t>/dx 
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This assumption, Eq. (^) is standard for simplifying the elliptic Helmhotz equation to an approximating 
parabolic equation. The resulting parabolic wave equation with reflection to first order in {(3 /a) is given by 

^ + i2fc[l - {2131a) cos(2fca;)] exp[i(2/3/a) sin(2A:a;)]|^ + {k^ - k'^)(j) = 



or 



^ +z2fcexp[i(2/3/a) exp(i2fca;)]^ + [k^ - k'^)<j) = 0. 
oz'' ox 
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Equation (0) is the parabolic equation that incorporates reflection. The difference between Eq. (|^) and the 
standard parabolic equation, Eq. (Q), is the additional factor exp[j(2/3/a) exp(i2fcx)] in the d(j)/dx term in 
Eq. (|^). Relative reflection as a function of the fraction (3 /a is thereby incorporated to first order into ^(a;, z). 

In order to account for the effect of reflection, contemporary solutions to the parabolic approximation 
solve the parabolkj-jsquation for an interacting pair of solutions (incident and reflected) or decouple the pair 
by simpliflcation.l30 Herein, we avoid the problem of coupled solutions. Our solution to Eq. (^ is a single 
synthetic wave that manifests both the incident and reflected wave throughout the domain. 

The initialization of (j) at some initial range, X;, over the depth column, z, renders the value, (f>{xi,z), 
over an open boundary thereby establishing the Dirichlet boundary conditions for a unique, stable solution 
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for 0J3 This initialization process is similar to that for the standard parabolic equation, but here we must 
also specify the fraction (3 /a. (As a starter, one could use UrickB and the references therein to predict (3 /a 
from reverberation and backscatter.) Solving Eq. (Q) (which is beyond the scope of this letter) in practice, 
one must not only take the usual precautions associated with the standard parabolic approximation but also 
take into account that Eq. d^) is an approximation that ignores some second-order terms of (P/a). 
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